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Abstract 

Background: Four human coronaviruses are currently known to infect the respiratory tract: 
human coronaviruses OC43 (HCoV-OC43) and 229E (HCoV-229E), SARS associated coronavirus 
(SARS-CoV) and the recently identified human coronavirus NL63 (HCoV-NL63). In this study we 
explored the incidence of HCoV-NL63 infection in children diagnosed with respiratory tract 
infections in Belgium. 


Methods: Samples from children hospitalized with respiratory diseases during the winter seasons 
of 2003 and 2004 were evaluated for the presence of HCoV-NL63 using a optimized 
pancoronavirus RT-PCR assay. 


Results: Seven HCoV-NL63 positive samples were identified, six were collected during January/ 
February 2003 and one at the end of February 2004. 


Conclusions: Our results support the notation that HCoV-NL63 can cause serious respiratory 
symptoms in children. Sequence analysis of the S gene showed that our isolates could be classified 
into two subtypes corresponding to the two prototype HCoV-NL63 sequences isolated in The 
Netherlands in 1988 and 2003, indicating that these two subtypes may currently be cocirculating. 


Background shaped spikes, giving the virus particle a characteristic 
Coronaviruses are large, enveloped, positive stranded §crown-like appearance. Three distinct groups of coronavi- 
RNA-viruses [1]. The viral RNA genome is 27-32 kb in _ ruses have been described based on serological affinity 
size, capped, polyadenylated and encapsidated inahelical and genome sequence. Coronaviruses can infect humans 
nucleocapsid. The envelope is studded with long, petal- and a variety of domestic animals and can cause highly 
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prevalent diseases such as respiratory, enteric, cardiovas- 
cular and neurologic disorders [2,3]. 


Until recently only three human coronaviruses were thor- 
oughly studied. Human coronavirus OC43 (HCoV-OC43; 
group 2) and human coronavirus 229E (HCoV-229E; 
sroup 1) were identified in the 1960s. They are responsi- 
ble for 10-30% of all common colds, and infections occur 
mainly during winter and early spring [4-7]. A third novel 
human coronavirus, SARS-CoV, was identified as the 
causal agent during the 2002-2003 outbreak of severe 
acute respiratory syndrome (SARS) [8-10]. Phylogenetic 
analysis showed that the SARS-CoV does not closely 
resemble any of the three previously known groups of 
coronaviruses, and therefore a tentative fourth group of 
coronaviruses was suggested [11,12]. However, an early 
split-off of the SARS-CoV from the coronavirus group 2 
lineage has also been suggested [13,14]. 


A new human coronavirus associated with respiratory ill- 
ness, HCoV-NLG3, was recently identified by a research 
team in The Netherlands [15]. The virus was isolated in 
January 2003 from a nasopharyngeal aspirate of a 7- 
month-old child suffering from bronchiolitis, conjuncti- 
vitis and fever. Screening of specimens from patients with 
respiratory symptoms identified seven additional HCoV- 
NL63 infected individuals, both children and adults, 
between December 2002 and February 2003. The com- 
plete viral gnome sequence was determined. The charac- 
teristic genome organisation of coronaviruses can be 
observed: the 5' two-third of the genome contains two 
large open reading frames (ORF), ORFla and ORF1Db. In 
the 3' part of the genome, genes encoding four structural 
proteins are found: spike (S), envelope (E), membrane 
(M), and nucleocapsid (N). The hemagglutinin-esterase 
(HE) gene, characteristic for group 2 coronaviruses, is not 
present in HCoV-NL63. Sequence analysis demonstrated 
that HCoV-NL63 shares 65% sequence identity with 
HCoV-229E. Phylogenetic analysis confirmed that HCoV- 
NLG3 is a new group 1 coronavirus, most closely related 
to HCoV-229E and porcine epidemic diarrhea virus 
(PEDV) [15]. Shortly after van der Hoek and colleagues 
published their discovery of the new human coronavirus 
HCoV-NL63, a second research group described the char- 
acterization of essentially the same virus [16]. The virus 
was isolated from a nose swab sample collected from an 
8-month-old child suffering from pneumonia in The 
Netherlands in April 1988. Real-time RT-PCR assays were 
designed for screening of respiratory tract samples. Four 
additional HCoV-NL63 positive samples, from children 
aged 3 months to 10 years, were detected between 
November 2000 and January 2001. 


HCoV-NL63 can be considered as a new important cause 
of respiratory illnesses and two different subtypes might 
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be currently cocirculating in the human population [15]. 
In this study we wanted to explore the incidence of HCoV- 
NL63 infection in children diagnosed with respiratory 
tract infections in Belgium. 


Methods 

Isolates and patients 

We studied 309 isolates from 279 patients with severe res- 
piratory symptoms collected from January 2003 until 
March 2004 at the University Hospital in Leuven, Bel- 
gium. These isolates originated from bronchoalveolar lav- 
ages, pharyngeal swabs, nasopharyngeal aspirates, and 
sputum samples. Routine diagnostic testing was per- 
formed for respiratory syncytial virus (RSV), influenza 
virus, parainfluenza virus and adenovirus. No _ prior 
amplification by cell culture was performed. The results of 
diagnostic tests for RSV were negative for 244 isolates, 
while 65 isolates were positive for RSV. Patients ranged in 
age from 1 month to 16 years, with a mean age of 2 years. 
The temporal distribution of the isolates corresponded to 
the yearly RSV epidemic period: 236 samples were col- 
lected from January to June 2003 and 73 samples were 
recovered during the first trimester of 2004 (Figure 1A). 


Pancoronavirus RT-PCR assay 

RNA was extracted from the collected specimens by using 
the QIAamp Viral RNA Mini kit (QIAGEN, Westburg, The 
Netherlands) according to instructions of the manufac- 
turer. Screening of the samples was performed by ampli- 
fying a 251 bp fragment of the polymerase gene using the 
following primer set: Cor-FW (5'-ACWCARHTVAAYYT- 
NAARTAYGC-3') and = Cor-RV_—(5'-TCRCAYTTDG- 
GRTARTCCCA-3') (Figure 2). These one-step RI-PCR 
assays (OneStep RT-PCR kit; QIAGEN) were undertaken 
in a 50 ul reaction volume containing 10 uL RNA-extract, 
10 ul 5x QIAGEN OneStep RT-PCR Buffer, 2 ul dNTP mix 
(final concentration of 400 uM of each dNTP), 1.8 pl QIA- 
GEN OneStep RT-PCR Enzyme Mix (a combination of 
Omniscript and Sensiscript reverse transcriptase and Hot- 
Starfaq DNA polymerase), 4 uM of each primer, and 
RNase-free water to 50 ul. The reaction was carried out 
with an initial reverse transcription step at 50°C for 30 
min, followed by PCR activation at 95°C for 15 min, 50 
cycles of amplification (30 sec at 94°C; 30 sec at 48°C; 1 
min at 72°C), and a final extension step at 72°C for 10 
min in a GeneAmp PCR system 9600 thermal cycler 
(Applied Biosystems, Foster City, CA, USA). PCR-products 
were run on a polyacrylamide gel, stained with ethidium 
bromide, and visualized under UV-light. 


RT-PCR assays for HCoV-NL63 

Samples that were found positive for HCoV-NLG63 were 
confirmed using one-step RT-PCR assays, which amplified 
four different regions of the HCoV-NL63 genome. Ampli- 
fication of a 314-bp gene fragment in the nucleocapsid 
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Detection of HCoV-NL63 and HCoV-OC43 in samples from patients suffering from severe respiratory symptoms. (A) 


Number of samples tested per month. (B) Patients infected with 
sample was isolated in April 2003 (not shown). 


HCoV-NL63 and HCoV-OC43. A single HCoV-229E positive 


region was performed with two specific HCoV-NL63 
primers: N5-PCR1 (5'-CTGTTACTTTGGCTTTAAAGAACT- 
TAGG-3', nt 26695-nt 26721) and N3-PCR1 (5'-CTCAC- 
TATCAAAGAATAACGCAGCCTG-3', nt 26982-nt 27008). 
Secondly a 237-bp fragment in ORF1b was amplified 


using the primers repSZ-1 and repSZ-3 described by van 
der Hoek and colleagues [15]. A third RT-PCR assay was 
carried out on the HCoV-NL63 positive samples amplify- 
ing a 839-bp fragment with ORFla specific primers: 
SS5852-5P and P4G1M-5-3P [15]. These one-step RT-PCR 
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Forward primer region 


HCoV-NL63 A CINCAIST GAATC 
HCoV-229E A CTR CAGHTINAATE 
FIPV A CHE CAAVATGAATER 
TGEYV ACE CAGHTAAATES 
PEDV NO AIAG CHICEWACC 
HCoV-O0C43 Lh CRCAMVATGAAREN 
BCoV Lh CRCAMVATGAAREN 
PHEV A CRCAMWAT GAA TER 
CRCV A CECAATGAATE 
MHV A CR CAINAT GAATS 
SDAV A CR CAINVAT GAAS 
SARS-CoV A CRCAINVAT GAAS 
IBV A CR CAMVATGAATEN 
TCOV A CRE CAIVAT GAAREN 
ACWCARHTVAAYYTNAARTAYGC 
Figure 2 
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Reverse primer region 


TGGGAYTAYCCHAARTGYGA 


Selection of primers for the novel pancoronavirus RT-PCR. Shown is the alignment of 14 coronaviral sequences of a conserved 
region of the polymerase gene. The forward (Cor-FW) and reverse (Cor-RV) primer sequences are shown at the bottom (Y = 
C/T, W = A/T, V = A/C/G, R = A/G, H = A/T/C, N = A/C/T/G). The coordinates of Cor-FW and Cor-RV are 14017 and 
14248, respectively, in the HCoV-NL63 complete genome sequence. The |4 coronavirus sequences used here are available 
from GenBank under the following accession numbers: HCoV-NL63, AY56/7487; HCoV-229E, AF304460; infectious bronchitis 
virus (IBV), Z30541; SARS-CoV, AY313906; HCoV-OC43, AY391777; PEDV, AF353511; bovine coronavirus (BCoV), 
AF391541; transmissible gastroenteritis virus, AF304460; MHV, X51939; PHEV, AFI 24988; sialodacryoadenitis virus (SDAV), 
AF 124990; turkey coronavirus (TCoV), AFI24991; canine respiratory coronavirus (CRCV), AY 150273; feline infectious perito- 


nitis virus (FIPV), AFI 24987. 


assays were performed essentially as described above. 
They were carried out using 5 uL RNA-extract and 0.6 uM 
of each primer. Only 45 cycles of amplification were run 
and annealing temperature was set at 50°C. Furthermore 
a 663 bp fragment of the spike gene was amplified using 
a RT-nested PCR. The outer primer set SINL5 (5'-GAGTTT- 
GATTAAGAGTGGTAGGTTG-3', nt 20391-nt 20415) and 
SINL3 (5'-AACAGTGTAGTTAACTACACGG-3', nt 21068- 
nt 21089) were used in a one-step RT-PCR, performed as 
described above, using 10 ul of RNA-extract and an 
annealing temperature of 48°C. A nested PCR was carried 
out with the inner primer set SINL5n (5'-GGTTGTTGT- 
TACGCAATAATGGTCGT-3', nt 20411-nt 20436) and 
SINL3n = (5'-ACACGGCCATTATGTGTGGTGAC-3', _ nt 
21051-nt 21073). The nested reaction mix was composed 
of 1 unit Tag polymerase, 1 ul of a 25 mM dNTP-mix, 10 
ul 5X buffer C (PCR Optimizer Kit, Invitrogen, The Neth- 


erlands), and 30 pmol of forward and reverse primer in a 
50 ul reaction volume. As template 10 ul of the outer PCR 
product was added. The cycling conditions were as fol- 
lows: an initial denaturation at 94°C for 5 min, followed 
by 40 cycles of amplification (45 sec at 94°C, 45 sec at 
54°C, 1 min at 72°C), and a final extension of 5 min at 
72°C. PCR-products were run on a polyacrylamide gel, 
stained with ethidium bromide, and visualized under UV- 
light. The amplicons were purified using the QIAquick 
PCR purification kit (QIAGEN) and sequenced with the 
respective primer pairs using the ABI PRISM BigDye Ter- 
minator Cycle Sequencing Reaction kit (version 3.1) on 
an ABI PRISM 3100 DNA sequencer (Applied Biosystems) 
according to the manufacturer's instructions. Positive and 
negative controls were included in each PCR experiment. 
The HCoV-NLG3 positive control was RNA isolated from 
a HCoV-NL6G3 culture. 
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Sequence analysis and phylogenetic analysis of the 
amplicons 

Chromatogram sequencing files were inspected with 
Chromas 2.2 (Technelysium Pty Ltd, Helensvale, Aus- 
tralia), and contigs were prepared using SeqMan II 
(DNASTAR, Madison, WI, USA). The obtained consensus 
sequences were compared with the prototype HCoV- 
NL63 sequences available in GenBank database release 
142.0 using BLAST analysis (NCBI BLAST server). Multi- 
ple sequence alignments were prepared using CLUSTAL X 
version 1.82 [26], and manually edited in the GeneDoc 
Alignment editor [27]. Phylogenetic analysis was con- 
ducted using MEGA version 2.1 [28]. 


Nucleotide sequence accession numbers 

The sequences determined in this study have been depos- 
ited in the GenBank sequence database under accession 
numbers AY758276 to AY758301. 


RT-PCR assays for HCoV-OC43 and HCoV-229E 

Our collection of samples was also screened using the 
pancoronavirus RT-PCR assay for the presence of HCoV- 
OC43 and HCoV-229E. Positive results were confirmed 
by one-step RT-PCR using HCoV-OC43 and HCoV-229E 
specific primer pairs located in the membrane glycopro- 
tein region (OC43-FW: 5'-GGCTTATGTGGCCCCTTACT- 
3', nt 28580-nt 28599; OC43-RV: 5'-GGCAAATCTI- 
GCCCAAGAATA-3', nt 28894-nt 28913; 229E-FW: 5'- 
TGGCCCCATTAAAAATGTGT-3', nt 24902-nt 24921; 
229E-RV: 5'-CCTGAACACCTGAAGCAAT-3', nt 25456-nt 
25475) [18]. One-step RT-PCR and sequence analysis 
were performed essentially as described above. Annealing 
conditions during the RT-PCR assay were modified: the 
annealing temperature was set at 55°C. 


Results 

Pancoronavirus RT-PCR assay 

A pancoronavirus RT-PCR assay is a usefull tool to test for 
all coronaviruses in a clinical sample. Besides quick 
screening for several pathogens in one assay, it supplies 
the possibility to identify previously unknown coronavi- 
ruses. The consensus RI-PCR assay as described by 
Stephensen et al., designed to amplify all known corona- 
viruses, is not able to detect HCoV-NL63 because of sev- 
eral mismatches with the primer sequences [15,17]. We 
modified these consensus primers based on an alignment 
of the HCoV-NL63 prototype sequence and 13 other 
coronavirus sequences (Figure 2). To determine whether 
the newly designed pancoronavirus RT-PCR assay effi- 
ciently amplifies a broad range of coronaviruses the RT- 
PCR assay was tested on cell culture supernatant of the 
four known human coronaviruses and three animal coro- 
naviruses: HCoV-NL63, HCoV-OC43, HCoV-229E, SARS- 
CoV, feline infectious peritonitis virus (FIPV), porcine 
hemagglutinating encephalomyelitis virus (PHEV), and 
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murine hepatitis virus (MHV). Amplification of the 
expected 251 bp region was observed for all tested corona- 
viruses (Figure 3). The sensitivity of the pancoronavirus 
RT-PCR assay was assessed by testing tenfold dilutions of 
HCoV-NL63 and HCoV-OC43 RNA. While 50 copies of 
HCoV-OC43 RNA copies per pl nasopharyngeal aspirate 
could be detected, the sensitivity for HCoV-NL63 was a bit 
lower i.e. 5 x 103 RNA copies per ul nasopharyngeal 
aspirate. 


Detection of HCoV-NL63 in clinical specimens 

The pancoronavirus RT-PCR assay was used for screening 
of specimens from hospitalized patients with respiratory 
symptoms collected between January 2003 and March 
2004. Samples, from which a 251 bp fragment could be 
amplified, were further identified by sequencing using the 
pancoronavirus primers. We studied 309 specimens with 
a temporal distribution that corresponded with the yearly 
RSV epidemic period (Figure 1A). A total of 244 samples 
were found negative for RSV by diagnostic testing. The 
279 patients in this study comprised of 211 patients aged 
<2 years (75.6%), 68 aged 2-16 years (24.4%). We 
detected HCoV-NLG3 in 7 samples (2.3%). One positive 
sample was collected at the end of January 2003 and coin- 
fection with RSV type B was present. Five of the positive 
samples were collected within a ten-day period at the end 
of February 2003, and one positive sample was collected 
at the end of February 2004, which showed coinfection 
with adenovirus and parainfluenza virus (Figure 1B, Table 
1). The seven positive samples were obtained from one 
patient aged 1 month, four patients of 1 year, one patient 
of 2 years, and one patient of 16 years. The patient files 
showed that all subjects suffered from respiratory tract ill- 
ness and some had underlying disease (Table 1). 


The seven HCoV-NL63 positive respiratory samples were 
confirmed by alternative RT-PCR assays. Amplification of 
a fragment of the nucleocapside gene and ORF1b was car- 
ried out. Sequence analysis of the N gene fragments and 
the ORF1b fragments showed 98-100% similarity to the 
prototype HCoV-NL63 sequences available in the Gen- 
Bank database (AY567487, AY518894). A third one-step 
RT-PCR was carried out for each positive sample to 
amplify part of the ORFla gene. Sequence analysis of the 
ORFla PCR-products revealed 99% sequence identity 
with both HCoV-NL63 prototype sequences available in 
GenBank. A neighbor-joining phylogenetic tree was 
constructed based on an alignment of the ORFla nucle- 
otide sequences from the HCoV-NLG63 positive samples 
and the available HCoV-NLG63 sequences in GenBank. 
HCoV-229E was used as an outgroup. The dendrogram 
shows that all HCoV-NL63 sequences cluster together, but 
two subclusters can be observed (Figure 4). 
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VI NL63 OC43 229E SARS PHEV FIPV MHV 


ss era 


Figure 3 
Gel electrophoresis after pancoronavirus RT-PCR assay. The indicated band of 251 bp corresponds with the expected ampli- 
con size. As a marker Molecular Weight Marker VI was used (Boehringer Mannheim, Germany). 


Table |: Patients hospitalized with respiratory tract illness associated with HCoV-NL63 infection 


Patient nr. Age Sex Symptoms Underlying disease Specimen Sample date 
| 153 | year male — URTI: fever, coughing, wheezing, rhinitis, diarrhoea none NPA 27 Jan 2003 
33545 I6 years male — LRTI: fever, coughing, respiratory distress, pharyngitis Smith-Lemli-Opitz syndrome NPA |4 Feb 2003 
21596 | year female LRTI: fever, coughing, respiratory distress Vater syndrome, epilepsy NPA 20 Feb 2003 
53887 | month female URTI: fever, rhinitis, two siblings have URTI none NPA 20 Feb 2003 
40001 | year male — LRTI: respiratory distress, cardiac arrest, rotavirus- epilepsy NPA 2| Feb 2003 
positive diarrhoea 

64880 2years male  URTI: fever, coughing, wheezing neurofibromatosis NPA 24 Feb 2003 
706885 | year female LRTI: pneumonia, fever, cyanosis, diarrhoea none PS 25 Feb 2004 


positive for RSV type B; »positive for adenovirus and parainfluenza virus 
LRTI, lower respiratory tract illness; URTI, upper respiratory tract illness; PS, pharyngeal swab; NPA, nasopharyngeal aspirate 
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Figure 4 
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Phylogenetic analysis of the partial ORFla nucleotide sequences. Accession numbers: HCoV-NL63, AY567487; HCoV-NL, 
AY518894; HCoV-229E, AF304460; NL-p466, AY567488; NL-p246, AY567489; NL-p251, AY567490; NL-p496, AY567491; 
NL-p223; AY567492; NL-p248, AY567493; NL-p72, AY567494; CAN39, AY675541; CAN52, AY675542; CAN57, AY675543; 
CAN 140, AY675544; CAN146, AY675545; CAN214, AY675546; CAN449, AY675547; CAN470, AY675548; CAN483, 
AY675549; CAN495, AY675550; CAN528, AY675551; CAN531, AY675552; CAN543, AY675553. 


Inspection of the two full gnome HCoV-NL63 sequences 
available in GenBank demonstrates that especially the 
aminoterminal region of the Spike protein can be very 
divergent. Therefore we decided to amplify this region to 
investigate the variability of these region in our patients. 
An RT-nested PCR assay was used to amplify part of the S 
gene. These partial spike sequences showed 98% similar- 
ity with the HCoV-NL63 prototype strains. An alignment 
of the S gene sequences from the Belgian samples, partial 
spike sequences from the positive samples identified in 
The Netherlands (data not shown), and the prototype 
HCoV-NL63 sequences, was used to constitute a neigh- 
bor-joining phylogenetic tree. The neighbor-joining tree 


was evaluated by 500 bootstrap pseudoreplicates. Two 
clusters can again be observed (Figure 5). 


Detection of HCoV-OC43 and HCoV-229E 

Screening of our sample collection for the presence of 
HCoV-OC43 and HCoV-229E was also performed. We 
detected HCoV-OC43 in 7 of 309 samples (2.3%) and 
HCoV-229E in one sample (0.3%). The seven HCoV- 
OC43 positive samples were collected during the winter 
and early spring of 2003 and 2004. The sample in which 
we detected HCoV-229E was collected in April 2003. The 
positive samples were confirmed by RT-PCR using specific 
HCoV-OC43 and HCoV-229E primer pairs that amplify 
part of the M gene [18]. The HCoV-OC43 and HCoV-229E 
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Figure 5 

Phylogenetic analysis of the partial S gene nucleotide 
sequences based on an alignment of the Belgian spike 
sequences, spike sequences from the positive samples identi- 
fied in The Netherlands, and the prototype HCoV-NL63 
sequences available in GenBank. Accession numbers: HCoV- 
NL63, AY567487; HCoV-NL, AY518894. PEDV was used as 
an outgroup. 


partial membrane sequences of the contemporary Belgian 
strains showed 97-99% similarity with the HCoV-OC43 
and HCoV-229E prototype sequences in GenBank. 


Discussion 

RSV, influenza viruses, adenoviruses, and parainfluenza- 
viruses are probably the most important viral agents of 
severe respiratory diseases. However, a substantial part of 
respiratory tract infections can not be attributed to any 
known pathogen. Underlying conditions and immuno- 
suppression largely determine the impact of respiratory 
viruses on individuals [19]. The common cold viruses 
HCoV-OC43 and HCoV-229E have also been associated 
with more severe lower respiratory tract conditions in 
infants and immunocompromised patients [20-23]. The 
clinical symptoms associated with HCoV-NL63 infections 
still need to be determined, but there are some indications 
that HCoV-NLG3 can cause severe respiratory illnesses in 
children and immunocompromised adults [15,16]. We 
detected HCoV-NL63, using a pancoronavirus RT-PCR, in 
patients suffering from relatively severe respiratory dis- 
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eases necessitating hospitalization. These positive sam- 
ples were collected from children aged 1 month to 16 
years. Two patients suffered from severe underlying 
disease: one patient suffered from Smith-Lemli-Opitz 
syndrome, a rare autosomal recessive disorder due to a 
primary enzymatic defect in the cholesterol metabolism. 
A second patient was diagnosed with VATER, a syndrome 
characterized by the sporadic association of specific birth 
defects or abnormalities such as vertebrae and vascular 
anomalies, anal atresia, trachea and esophagus problems 
and renal anomalities. All HCoV-NLG3 infected patients 
established a complete recovery from their respiratory 
symptoms. One-step RT-PCR assays were used to detect 
and confirm these positive samples. 


Results from epidemiological surveys conducted in the 
1970's have led to the conclusion that human coronavi- 
ruses are distributed worldwide and circulate during sea- 
sonal outbreaks [22]. Our results indicate that HCoV- 
NLG3 is the causal agent in a significant portion of respi- 
ratory diseases of unknown etiology. We detected HCoV- 
NL63 in respiratory samples collected in February 2003, 
with a frequency of 7.1%, and during February 2004, with 
a frequency of 2.5%. These results seem to support the 
tendency of human coronaviruses to circulate mainly dur- 
ing the winter season |7,24]. However, in this study, 
sampling was only performed from January to May during 
the yearly RSV epidemic period, while no samples from 
the summer and autumn months were screened. The first 
publication on HCoV-NL63 showed that the virus circu- 
lated in Amsterdam during the winter months of 2002/ 
2003 [15]. More recently another set of Amsterdam sam- 
ples was screened, obtained during the winter of 2001/ 
2002 and 2003/2004. HCoV-NL63 was found in one tra- 
chea sample obtained in February 2002, and in two 
oropharyngeal aspirates from December 2003 and Janu- 
ary 2004, respectively (data not shown). Combined with 
the data that we present here from Belgium, these findings 
confirm that HCoV-NL63 reappears each winter season 
similar to the previously known respiratory viruses. 
Recently, research teams from Australia, Japan and Can- 
ada, have submitted partial HCoV-NLG63 sequences to the 
GenBank database (AY600442-AY600446, AY662694- 
AY662698, AY675541-AY675553). This indicates that 
this newly discovered human coronavirus has a world- 
wide distribution. 


Sequence analysis of the highly conserved nucleocapsid 
region showed that the Belgian isolates are similar to the 
two prototype HCoV-NL63 complete genome sequences 
in GenBank isolated in the Netherlands in 1988 and 
2003. Furthermore, phylogenetic analysis of part of the 
ORF 1a region of our patients showed the same subclus- 
ters of HCoV-NL63 that were described previously [15] 
(Figure 4). This finding supports the suggestion that 
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several HCoV-NL63 subtypes with distinct molecular 
markers are cocirculating, also in Belgium. A large insert 
in the 5' part of the S gene of HCoV-NL63 compared with 
HCoV-229E has been described [15,16]. Both HCoV- 
NL63 complete genome sequences show only 89% 
sequence identity in this spike insert region, which 
implies that there are at least two different HCoV-NL63 
subtypes. Sequence analysis of this spike insert region 
revealed that our samples show similarity to both proto- 
type HCoV-NLG63 subtypes, which was confirmed by phy- 
logenetic analysis. The partial S gene sequences cluster 
together with the two prototype HCoV-NL63 sequences in 
two different groups (Figure 5). This confirms that the 
HCoV-NLG3 subtypes first isolated in 1988 and 2003 are 
cocirculating. When analysing the dendrograms based on 
ORF1a and S gene sequences, a discordance in the cluster- 
ing pattern of some HCoV-NL63 isolates (e.g. HCoV-NL 
and NL-p223) can be observed, suggesting a possible 
recombination event. Further research of complete 
genome sequences of these isolates is required. Drawing 
conclusions based on phylogenetic analysis of one single 
gene therefore requires caution as the true phylogeny can 
only be demonstrated by analysing complete genome 
sequences. 


Screening of our sample collection for the presence of 
HCoV-OC43 and HCoV-229E revealed seven HCoV- 
OC43 positive samples and only one HCoV-229E positive 
sample. All positive samples were isolated during winter 
and early spring, which is concordant with the results of 
previous epidemiological studies. HCoV-OC43 infected 
samples were mainly identified during February 2003 and 
February 2004 (Figure 1B). These data show that the epi- 
demic seasons of HCoV-OC43 and HCoV-NL63 coincide. 
The positive samples were collected from children aged 1 
to 12 years, whom all suffered from respiratory symp- 
toms. The very low detection rate of HCoV-229E com- 
pared with the frequent detection of HCoV-NL63, might 
imply that HCoV-NL63, closely related to HCoV-229E, is 
currently more important as a causal agent of respiratory 
diseases. At the moment, there are no data concerning 
cross-neutralization between HCoV-229E and HCoV- 
NL63. In theory, such cross-neutralization might be pos- 
sible, since both viruses are relatively closely related spe- 
cies belonging to coronavirus group 1. Antigenic cross- 
reactivity has already been demonstrated between SARS- 
CoV and group 1 coronaviruses TGEV, FIPV and CCoV 
[25]. 


The development of a pancoronavirus RT-PCR assay using 
a primer set that matches all known coronaviruses might 
be useful for the identification of new coronaviruses. This 
pancoronavirus RT-PCR-assay can also be used as a diag- 
nostic tool to detect any of the four currently known 
human coronaviruses in clinical samples. 


http:/www.biomedcentral.com/1471-2334/5/6 


Conclusions 

Human coronavirus NLG3 is a new important respiratory 
pathogen that can cause severe respiratory infections in 
children. Sequence analysis of the HCoV-NL63 isolates 
detected in our study demonstrates that our Belgian iso- 
lates can be classified into two subtypes corresponding to 
the two prototype HCoV-NL63 sequences isolated in The 
Netherlands in 1988 and 2003. Our findings indicate that 
these two subtypes may currently be cocirculating. 
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